Modelling a suspended nanotube oscillator 
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We present a general study of oscillations in suspended one-dimensional elastic systems clamped 
at each end, exploring a wide range of slack (excess length) and downward external forces. Our 
results apply directly to recent experiments in nanotube and silicon nanowire oscillators. We find 
the behavior to simplify in three well-defined regimes which we present in a dimensionless phase 
diagram. The frequencies of vibration of such systems are found to be extremely sensitive to slack. 



Vibrations of one- dimensional systems (i.e. systems 
with cross-sectional dimension much smaller than their 
length) suspended under the influence of a downward 
force have long been of interest in the context of such ap- 
plications as beams, cables supporting suspension bridges 
and ship moorings 0, S 01- Such systems display a 
wide range of behavior, depending on the amount of 
slack present in the system, the downward force, and 
the aspect ratio. Previous studies, which have largely 
been analytical, have been restricted to certain limit- 
ing cases of these parameters. Recently work on oscil- 
lating nanoscale systems, such as carbon nanotubesQ 
and silicon nanowires^, has opened the possibility of 
experimentally exploring such vibrations in entirely new 
regimes. In the present work, we study numerically the 
oscillations of a one dimensional elastic system over an 
extensive range of both the slack and force parameters, 
providing insight into the physics which separates this 
parameter space into three distinct regimes of behavior. 

The treatment below is entirely general with illustra- 
tive examples taken from parameters relevant to carbon 
nanotubes. Since their discovery in 1991 3, nanotubes 
have found many applications in device technology due 
to their small size, robust structure and superior elastic 
properties [lllalSlUil- Many of these applications involve 
the use of nanotubes as mechanical oscillators, making 
theoretical understanding of the vibrational properties of 
nanotubes in various geometries of current interest [3. 1 lll|. 
Recent experiments |l2l | have studied the behavior of the 
transverse vibrations of a suspended nanotube clamped 
at both ends, under the action of a downward force, as 
sketched in Figure ^ These suspended nanotubes gen- 
erally have around 1% slack, denoted in this work by s, 
which we define to be the ratio of the excess length of 
the tube to the distance between clamping points. 

Analytic model and computational techniques — The 
potential energy of a one-dimensional elastic continuum 
under a uniform downward force is 




Figure 1. (a) Suspended nanotube under the action 
of a gate voltage V. (b) Jump-rope mode 



respectively; E and F are the extensional and flexural 
rigidities; and / is the downward force per unit length. 
For single-walled nanotubes with diameters d > I nm, 
the rigidities we obtain using the Tcrsoff-Brenner empir- 
ical potential [l3| are, to within 0.5%, 
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regardless of the chirality of the tube [l^. We note that 
essentially the same values of E and F can be obtained 
using the in-plane elastic modulus of graphite. 

To evaluate the potential energy (Equation^, we be- 
gin with r(Z), the vector location of the segment at posi- 
tion I from the end of the tube. The local strain is, then. 
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u{l) ^ |r'| - 1 
and inverse radius of curvature is 



Here, I represents distance along the system of un- 
stretched length L; u{l), R{1) and z{l) represent the lo- 
cal strain, radius of curvature and vertical displacement, 
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where prime indicates differentiation with respect to /, 
and the tangent vector i = r'/ |r'|. 

For a given force and slack, we discretize r(^) along I 
with a number of points (typicaUy 30) sufhcient to render 
all plots shown below visually indistinguishable under in- 
creasing resolution. To find the equilibrium configuration 
of the system under the action of the external force, we 
then relax the system using a conjugate gradients algo- 
rithm with forces computed using finite differences. Fi- 
nally, we compute the force constant matrix Kij = q^ q^ 
and diagonalize it to obtain frequencies v while assigning 
a mass /i to each discretized point, where /i and AZ are 
the mass per unit length of the system and the distance 
between sample points respectively. 

For a one-dimensional system, slack s, length L and 
downward force per unit length / can be regarded as ex- 
ternal parameters that can be changed independently. 
Mass density /i and rigidities E and F on the other 
hand, all depend on the cross-sectional geometry of the 
system and can therefore be considered internal param- 
eters. To make our discussion applicable to a general 
one-dimensional system, we present our results in units 
scaled by appropriate combinations of these basic param- 
eters. In addition, we give results in physical units for 
the concrete system of a nanotube of diameter 2 nm with 
length L = 1.75 /xm. For frequency plots, besides the 
physical frequency we plot a dimensionless frequency 
D = v^{lxL^/F). 

As a dimensionless measure of this force, we choose a 
force control parameter V = \J fL'^/F, proportional to 
the square root of the force. In recent experiments on 
suspended nanotubes^3|j the external force was intro- 
duced via a gate voltage, as shown in Figure Such 
a voltage results in a force per unit length / = \^V'^ . 
Here c is the capacitance per unit length of the nanotube 
given by 
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where eo is the dielectric constant of vacuum and z is 
the distance between the wire and gate. For the concrete 
example of a L=1.75/xm nanotube with a diameter of 
2nm, we take the tube to be suspended 500nm above the 
gate. In this dimensionless parameter of ^ = 10 

corresponts to a voltage of . As another example, 
a 100/im-long Si nanowire of lOOnm diameter, with a z 
of 25^m would require a voltage of ~60V to attain the 
same force control parameter V of ten. Moving to larger 
systems, a 0.6mm-diameter steel piano wire which is ^2m 
long under the influence of gravity would have the same 
dimensionless parameter of V = 10. 

Figures |H and |31 display our primary results, both in 
dimensionless form and for the nanotube example, for 
the dependence of the vibrational frequency on slack 
and force. Figure |21 shows the dependence of frequency 
on slack over the range of small slacks for zero force. 
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Figure 2. Frequency versus slack at zero force control 
parameter: classic solution for zero slack (points on 
vertical axis), numerical results for non-zero slack 
(curves). The dimensionless frequency v — v^JJlLAJF 

applies to a general system. The modes are 
doubly-degenerate and odd (sohd curves), out-of-plane 
and even (dashed curves) or in-plane and even (dotted 
curves). Arrows indicate mode crossings. 



The style of the curves indicate whether the mode is 
doubly-degenerate and odd (solid), out-of-plane and even 
(dashed) or in-plane and even (dotted). After the mode 
crossings shown with arrows in Figure 13 all curves re- 
main flat up through a slack of 2%, the largest slack 
considered in this paper. Figure 01 shows the dependence 
of the frequency on the force control parameter. Fig- 
ure|2ta) shows the symmetry of the modes as in Figure|21 
for a system of 1% slack. Figure Olb) shows slacks of 
5=0.25%, 1% and 2% with the horizontal axis scaled by 
^ which, for reasons described below, causes the curves 
to collapse upon each other for all but the smallest slack. 

We find three simplifying limits for the physical behav- 
ior of the system as a function of slack s and force control 
parameter v. In each such limit, a different term in the 
potential energy (Equation ^ dominates with the other 
terms either being negligible or acting as a constraint. 
These three regimes, which we term buckled beam, hang- 
ing chain and hanging spring, are shown in Figure 01 

For small external forces, the extensional rigidity acts 
as a constraint which keeps the tube length fixed. If 
the force is sufficiently small, the bending term in Equa- 
tion dominates the dynamics and the system acts like 
an Euler-buckled beam. The in-plane vibrations of buck- 
led beams have been studied analytically, but only in the 
absence of an applied force and only under conditions of 
very smaU slack (under - 0.07% 0). Our results extend 
this treatment to three dimensions and larger values of 
slack, probing at the same time the dependence on the 
downward force. In the absense of applied force all modes 
are doubly degenerate at zero slack with the frequencies 
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Figure 4. The three regimes in which nanotube 
behavior simplifies, as a function of voltage and slack. 
The dashed line between the "hanging spring" and 
"hanging chain" regimes is the crossing between the 
second even in-plane mode and the second odd in-plane 
mode. The dotted line separating the "buckled beam" 
regime from the "hanging chain" regime is the contour 
on which the derivative li2Si£. j^as a value of —1/8. The 

a log s ' 

axis on the right shows the force control parameter, 



overlaying the slack axis of Figure El No such symmetry 
applies to other out-of-plane modes, whose frequencies 
drop slightly, and then remain essentially independent of 
slack. 



Figure 3. (a) Dimensionless frequency versus force 
control parameter: in-plane modes (dashed curves), 
out-of-plane modes (dotted curves), degeneracy of in- 

and out-of- plane modes (sohd curves), (b) 
Dimensionless frequency versus force control parameter 
scaled by fourth-root of slack {v/ ^s/1%) for three 
values of slack s: 0.25% ('+'), 1% (curves), 2% ('©'). 

Frequency in MHz and voltage in V refer to the 
nanotube example described in the text. The scaled 

frequency iy ~ v^JflLA/F and the force control 
parameter v — ^ fl?jF apply to a general system. 

of the classic unbuckled clamped beam shown as 
circles on the vertical axis of Figure HJ 

With the introduction of slack, the lowest out-of-plane 
mode, which we term "jump rope" , corresponds to rota- 
tion of the relaxed system about the clamping points, as 
illustrated in Figure Rotational symmetry requires 
that this mode has zero frequency in the absence of an 
external force. It therefore shows up as a dashed line 



In contrast to the out-of-plane modes, the in-plane 
modes are affected by a length-conservation constraint 
which changes with the amount of slack. For reasonably 
small slack, the odd modes (shown as solid lines) are un- 
affected by this constraint, since they conserve length by 
symmetry, and thus remain degenerate with their out-of- 
plane counterparts. The even modes (shown as dotted 
lines), in order to conserve the tube length, acquire two 
extra nodes as slack is introduced, leading to a mode 
crossing between the even and odd modes (shown with 
arrows in Figure EJ. 

When the force is large enough to overcome the bend- 
ing term in Equation but still too small to stretch the 
nanotube, the tube behaves as a hanging chain clamped 
at both ends, whose equilibrium profile forms the well- 
known catenary. In this regime, the slack dependence of 
the frequency can be solved for analytically in the limit 
of small slack. In the small slack limit, the shape of the 
tube is parabolic, and the tension is constant throughout 
the tube. The frequency is then determined by the result 
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for a simple string under tension, giving 

uin — nn. (6) 

Since in this equation the frequency is proportional to 
s'-^^^, we scale the horizontal axis in FigureOIb) by y^, 
causing the curves for each slack value to collapse for 
small forces, although for 0.25% slack there is deviation 
in the higher modes, in particular in the mode pointed 
to by the arrow. 

For a given force, as the slack decreases, the tension 
increases and eventually the tube begins to stretch and 
departs from the hanging chain regime. This avoids the 
infinite frequency at very small slack which is predicted 
by Equation El This effect is responsible for the devia- 
tion from linear behavior of the tube with 0.25% slack in 
Figure 13 in contrast to the tubes with 1% and 2% slack. 

In the large force, small slack regime, the extensional 
term in Equation ^ can no longer be treated as a con- 
straint, as the tube begins to stretch. The nanotube 
then behaves as a hanging spring clamped at both ends. 
The frequency in this regime can be derived in a manner 
similar to that in the hanging chain regime: the shape 
of the nanotube is parabolic and the tension is constant 
throughout the tube. However, in the case of the hang- 
ing spring, the tension is proportional to the change in 
length of the tube. Under these conditions, the frequency 
is proportional to the cube root of the force: 

The frequency is independent of slack, because as the 
tube is stretched the original amount of slack in the nan- 
otube becomes insignificant. However, the point of onset 
of the hanging spring regime does depend on slack. 

In the buckled beam and hanging chain regimes, the 
length conservation constraint in a slack tube forces the 
first even in-plane mode to have two nodes. In contrast, 
in the hanging spring regime there is no length conserva- 
tion constraint, and consequently, the lowest frequency 
in-plane mode has no nodes. As the transition between 
these regimes occurs, there is a mode crossing between 
the even and odd in-plane modes, shown with an arrow 
in Figure 131 in the curve for 0.25% slack. The crossing 
between the second even in-plane mode and the second 
odd in-plane mode is shown as the dashed curve in Fig- 
ure01which separates the hanging spring regime from the 
hanging chain regime. 

In conclusion, we have studied numerically and ana- 
lytically, as a function of both the amount of slack and 
applied force, the transverse vibration frequencies of one 
dimensional suspended systems clamped at both ends. 



We have chosen as a concrete application to focus on sus- 
pended nanotubes under the influence of a gate voltage 
which have recently been studied experimentally ■ We 
find three regimes in which the behavior simplifies and 
can be understood using approximate analytic solutions. 
This understanding is important in characterizing and 
designing nanoresonators where it can prove difficult to 
control or measure parameters such as slack. The hang- 
ing chain regime is particularly suitable for tuning and 
characterizing nano-oscillators, since Equation El which 
governs the frequency in this regime, remains valid inde- 
pendently of imperfections affecting the boundary condi- 
tions at the clamps. 

This work was supported by the NSF through the Cor- 
nell Center for Materials Research and the NIRT pro- 
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